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Abstract 

The so-called Mr scaling for dilepton spectra at RHIC energies is ex- 
amined. It is seen that a proper accounting of the complete set of dilepton 
producing processes forces us to abandon the proposed scaling around Mr = 
2.6 GeV, which has been put forward as a possible signature of the presence 
of quark-matter. Any substantial transverse expansion in the QGP phase 
itself will also offset this scaling behaviour. The rates of lepton pair produc- 
tion are time-integrated in Bjorken hydrodynamics and the different sources 
are compared against each other. 
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1 INTRODUCTION 



Quantum Chromodynamics is by now firmly established as the theory of 
strong interaction. One of its most spectacular predictions, namely the for- 
mation of quark-gluon plasma (QGP) will be experimentally investigated at 
the Relativistic Heavy Ion Collider under construction at Brookhaven and 
the proposed Large Hadron Collider at CERN. The QGP likely to be pro- 
duced in these experiments is expected to survive only for a brief duration of 
at the most a few fm/c, and a host of suggestions have been put forward to 
verify its formation and subsequent evolution. The early theoretical descrip- 
tions |I| have been tremendously successful in generating activity and ideas. 
Recently, approaches have been put forward to treat certain aspects of the 
nucleus-nucleus collisions on a microscopic level 0. If the QCD plasma is 
formed in such collisions and if QCD as a theory admits a first order phase 
transition, the system might further pass through a phase mixture of quarks, 
gluons and hadrons 0. Subsequently, the hadrons lose thermal contact and 
free-stream towards the detectors. 

The space-time evolution of this matter could be reflected by the yield 
of high energy photons and dileptons. These electromagnetic probes are 
considered reliable because their production rates are very strongly increasing 
functions of the temperature (local energy density) and they escape from the 
system without attenuation after they are produced. 

These endearing aspects of photons and dileptons have led to a thorough 
examination of the processes affecting their production. Thus, for example, a 
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study of photon-yielding processes in hadrons and QGP led to the surprising 
finding that the rate of emission of photons at a given temperature is phase 
independent This necessitated much higher initial temperatures for the 
QGP phase for cleaner signals to emerge 0. Such high initial temperatures 
may perhaps be achieved at the RHIC or the LHC ||. 

For a long time the dileptons having large invariant masses were believed 
to have their origin in the process qq — > 11 for the plasma and in the process 
7r + 7r~ — > 11 for the hadronic matter. As the pion-annihilation process thus 
envisaged is dominated by the p meson, it immediately led to the suggestion 
that dileptons having their invariant masses larger than that for the p meson 
would have their origin predominantly in the quark matter ||. 

These early expectations need to be re-examined in the light of recent 
findings |7] , where for the first time a complete list of light mesons was used 
to obtain the rates for dilepton yield in a hot hadronic matter. An important 
consequence of the quoted work is that the inclusion of vector decays and two- 
body reactions lead to an enhancement of the rate for dilepton production by 
a factor of 10-50, compared to the situation when only pion annihilation is 
considered, up to an invariant mass M ~ 3 GeV. An immediate outcome of 
this result is that the safe window for detecting dileptons from the plasma is 
now pushed up to a region where Drell-Yan dileptons are no longer negligible. 

These results also force us to abandon the so-called scaling for dilep- 
ton spectra around = 2.6 GeV, which has been proposed gnature 
of the presence of quark-matter [|J. The transverse mass scaling of the lepton 
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spectrum in high energy heavy ion collisions will occur for an ensemble of 
massless quarks in thermal equilibrium, undergoing longitudinal expansion 
only. These conditions ensure that the only scale in the problem is the tem- 
perature, T. Then dN / dM 2 dyd 2 qT is a function of Mt only. However, this 
feature does not persist in the hadronic sector where the timelike electro- 
magnetic form factors explicitly break the scaling law. Nevertheless, there 
has been recent hope that the scaling in the pure QCD sector could still 
be observed in the net signal ||. As we shall see, this possibility hinges on 
the dilepton yielding processes in the nonperturbative regime being approx- 
imated by 7r + 7r~ — > 11, which has a negligible contribution beyond M > 1 
GeV. Unfortunately, the Mt scaling is not observed even for a purely longi- 
tudinal expansion of the system once a full set of dilepton producing rates in 
hadronic matter are accounted for. Also, we mention the unavoidable com- 
plications resulting from transverse expansion ||. Even in the absence of the 
arguments made in this paper, this would seriously challenge the extraction 
of a clean signal. 

In this paper, we consider the impact of a larger set of lepton pair emitting 
processes on lepton spectra and on the hypothesis of Mj> scaling. 

2 FORMULATION 

The dilepton yield calculations for collisions involving two nuclei have rou- 
tinely been carried on in the following manner. First, the Drell-Yan con- 
tribution is either calculated or estimated by scaling the pp spectrum by 



4 



an appropriate power of A. A second contribution is obtained by integrat- 
ing the thermal emission rate over the space-time history of the collision 
|J. We assume that a thermalized quark-gluon plasma is formed at some 
initial temperature To at initial proper time tq. This plasma is then be- 
lieved to expand and cool. Soon it reaches the critical temperature T c for an 
assumed first-order chiral symmetry /deconfinement phase transition. The 
system is assumed to continue to expand at this fixed temperature, deriving 
the necessary energy from the latent heat of the transition. Eventually the 
quark-gluon plasma is entirely converted to hadronic matter. After this, the 
hadrons might maintain thermal contact for a while, as the matter continues 
to expand and cool. Finally, one estimates a 'freeze-out' temperature Ty at 
which the hadrons lose thermal contact with each other. Then they begin 
free-streaming towards the detectors. One adds to the Drell-Yan contribu- 
tion the emission over the cooling curve, from To to Tf, using equations of 
state for quark-gluon plasma and hadronic gas. This is the way in which the 



dilepton mass spectrum was computed in ref. [10|, for example. 



2.1 The rates for dilepton production from the QGP 
and hot mesons 

Dileptons are normally assumed to be produced via the reaction qq — ► 
in the plasma phase and 7r + 7r~ — > in the hadron phase. In the QCD 
sector we shall ignore corrections to the Born term. Those effects, and their 
consequence on Mt scaling are worthy of interest. Thus, at this level of 
approximation, the corresponding rate for production of dileptons having an 
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invariant mass M from the QGP phase is known to be given by 



dN a q {M) 
d 4 xdM 2 ~ 2 (2tt) 4 



M 3 TKi{M/T) 



Am 2 



(1) 



where m q is the mass of the quarks undergoing annihilation, 



a q {M) 



An a 2 

YJI 2 



1 + 



2mf 

H 2 
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1/2 



F g (M) 



and, 
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F q (M) = N c (2s + 1) £ 4 - - 



(2) 



(3) 



for the QGP consisting of u, d, and s quarks. In the above, m; is the mass 
of the leptons, a is the electromagnetic fine structure constant and N c is a 
color factor. The quark and lepton masses will ultimately be set to zero. 

The corresponding rate for pionic annihilation is obtained by replacing 
m q by m n in Eq. ([I]) and further taking the pionic form-factor F n in place 
of F q in Eq. For the pion electromagnetic form factor we use a recent 



parametrization of Biagini et al. [11 



We have already indicated that the rates of dilepton yielding processes 
in a hot meson gas have recently been computed [0 by taking a complete 
set of light mesons and including decays and two-body reactions. Because of 
phase space arguments, these contributions should be dominant. The light 
pseudoscalar (P) and vector (V) mesons we included consist of ir, r], p, u, 
1]', <f), K, and K*. Below the na\ — > e + e~ threshold, we have considered 
the decays p ne + e~, K*± -> A^e+e - , K*°(K*°) -> K°(K°)e + e-, 

<■ r]e + e~, r]' — > p°e + e~, r( — > ue + e~, (ft r/e + e~, 



uj — > 7r°e + e , p° 
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— > r]'e + e~ , — > 7r°e + e~. We include all possible reactions of the type P 
+ P — > e + e~ and V + V — > e + e~. The V + P initial states considered are 
u 7T°, p 7r, (f) 7T°, u) i], <j) rj, p° T], to T)' ', 7/', p t/', and K*K. Furthermore, 
the effect of the ai meson has been discussed in connection with real photon 
emission 



121, |13|1 . We have added the rates corresponding to the initial states 



nai and a\a\ into electron-positron pairs. The formalism to calculate all 
the Feynman amplitudes above is that of effective chiral Lagrangian models 
[0, 14]. Ref. @ contains a complete analysis of the relative contribution of the 
processes we have included. The reader will also find there a discussion of how 
our Lagrangian parameters have been adjusted to reproduce experimental 
data and of form factor issues. 

We have found that the net rate obtained by summing all the processes 
scales as TK\(M/T) to a high degree of accuracy between M=0.3 GeV, and 
3 GeV. This is because of the fact that the majority of reactions considered 
in Ref. @ share the pion electromagnetic form factor. In order to see this 
most clearly, we have plotted in Fig. 1 the effective form-factor obtained by 
summing all the channels enumerated above as, 



dN 



d A xdM 2 2 (2tt 



aeff(M) M 3 TK!(M/T) 



4m; 



with 



47T a 2 



1 + 



2mf 



1 - 



4mf 



1/2 



Peff(M) 



(4) 



(5) 



at temperatures of 100, 150, and 200 MeV. The scaling with TK^M/T) 
mentioned above is seen to be satisfied over a reasonable range of invariant 
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mass. The deviation seen at very low-masses has its origin in the contribu- 
tions lo —>■ tc°11, p — > nil, and <fi — > tt°11, which however are only marginal for 
large invariant masses. There is a small breaking of our scaling hypothesis 
above M = 2.5 GeV because of the a%a\ — > e + e~ threshold. We have also 
shown for comparison the form-factors for the pionic annihilation and the an- 
nihilation of quarks. It is seen that if only pionic annihilation processes are 
included, the quark annihilation processes will dominate the spectra beyond 
M pa 1.5 GeV, by two orders of magnitude. However, a proper accounting 
of all the mesons present in the hot hadronic matter makes the identification 
of the quark-gluon signal considerably more difficult; we shall return to this 
point later. 

2.2 The dilepton distributions 

2.2.1 Results for transverse-mass distributions 

The simplest space-time evolution dynamics is provided by the longitudinal 
hydrodynamic expansion model of Bjorken 0. However after the proper 
time ph R/c s , where R is the transverse dimension of the system and c s is 
the speed of sound, the rarefaction wave-front from the surface will reach the 
centre and the transverse expansion |§ of the system can not be ignored. The 
transverse velocity of the fluid elements would increase monotonically if there 
is no change in equation of state. However, a first order phase transition, as 
in the case of quark-gluon plasma, slows down the transverse expansion of 
the system considerably. 

We have considered head-on collision of two gold nuclei, and used the 
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criterion 

, 2vr 4 dN 

~ 45C(3)7r J R 2 4a Q dy n 

to relate initial temperature (T ) and the time (tq) to the pion rapidity den- 
sity (dN/dy^). In the above, /o=0 gives the fraction of the QGP at the initial 
time, and we have taken oq = 47.57r 2 /90 for a system consisting of u, d, and 
s quarks. The fraction f = when T < T c , < /o < 1 when T = T c , and 
/o = 1 when Tj > T c . Furthermore, r, the ratio of degrees of freedom between 
the QGP phase and the hadronic phases is 47.5/6.8, which corresponds to 
hadronic matter having the complete set of light mesons tt, p, uj, tj, rj' , 0, K, 
K*, and a\ 0. As the immediate purpose of this work is to examine the rel- 
evance of Mt scaling, we shall use Bjorken's hydrodynamic expansion model 
in the following. Now the volume element in Bjorken hydrodynamics is given 
by 

d 4 x = d 2 %T dz dt = iiR 2 dy rdr (7) 

In order to see the so-called scaling behaviour of the Mt spectra, we first 
consider with To =250 MeV and tq=1 fm/c. With those values, the 

system is initially in the QGP phase during the proper time from tq to 
r q = (T /T c ) 3 r , in the mixed-phase during r q to = rr q , and in the hadronic 
phase during to 77 = (T c /Tf) 3 Th, beyond which it undergoes freeze-out. 
We have taken the critical and freeze-out temperatures as 160 MeV and 120 
MeV, respectively. 

Now the number of dilepton pairs having an invariant mass M and trans- 



verse momentum p?, so that Mt = \/Pt + M 2 , can be written as (see [ITU 
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for comparison), 

dN 
dM 2 dM T dy 



a q (M)M 2 M 7 
4(2tt) 4 
a q (M)M 2 M T 



4(2tt) 4 
a cS (M)M 2 M T 

4(2tt) 4 
a cS (M)M 2 M T 

4(2vr) 4 



1 _ 
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i _ 4m 2 
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M 2 
4m 2 
M 2 " 
4m 2 
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r - 1)t 2 ttR 2 K (M t /T c ) 



-r{r - 1)t 2 ttR 2 K (M t /T c 
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R 2 3T^t^Mt 6 [G(M t /T c ) - G{M T /Tf)\ 



where 



G(z) = z 3 (8 + z 2 )K 3 (z) . 



(9) 



In Fig. 2 we have plotted the results for the transverse mass distribu- 
tion of dileptons for a fixed My = 2.6 GeV ||, against p?. Note that this 
way of plotting is convenient: if the scaling holds this plot will be a straight 
horizontal line. We show results obtained when only the pionic annihilation 
process is included for the hadronic matter, and when the entire spectrum 
of hadronic processes leading to dileptons are considered. If we restrict our 
analysis to the pion-pion process, the total contribution obtained by sum- 
ming up the quark matter and hadronic matter signals will exhibit scaling. 
However, taking account of all the processes in the hadronic matter leads to 
a considerable enhancement of dilepton production even up to an invariant 
mass of 2 GeV (lower values of px on this figure) f7j. No scaling of the M T 
spectra will be seen as the QGP contribution, which itself scales with 
for the case of no transverse expansion, is only a small fraction of the total 
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contribution. In this sense, the apparent flatness below pt = 1.5 GeV/c is 
accidental and not a reflection of the QCD M T scaling. The bulk of the signal 
is from the hadronic sector, where the hadronic electromagnetic form factors 
break the scaling. It should also be remembered that, even if the lifetime of 
the QGP is large, as it would be in the case of a higher initial temperature, 
the transverse velocity of the fluid could then be substantial even during the 
QGP phase which would furthermore offset the M T scaling property even for 
the QGP contribution. 

Thus we conclude that the M T scaling of the dilepton spectra, which has 
been proposed as a signature of the presence of quark-matter, does not hold at 
the temperatures considered in this work, once a more complete description 
of dilepton-yielding processes is achieved. We do not display in this work 
the Mt feature of the Drell-Yan dileptons. We shall return to this point in 
the discussion. 

2.2.2 The invariant mass distribution 

The invariant mass distribution of dileptons provides a very useful parametriza- 
tion of the data, as they are relatively less affected by the models for space- 
time evolution. The final results are easily written for the initial state in 



QGP as, 



dM 2 dy 



dN 



2(2tt) 4 




<J q (M) 
2(2tt) 4 
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Te<E(M) M3 



2(2vr) 



2(2tt) 



1 - 
1 



4m; 
Am 2 

TV 



irR 2 -r{r -1)t 2 T c K x {M/T c ) 
ttR 2 3T^t^M- 5 [H(M/T c ) - H(M/Tf)] 



(10) 



where 

H{z) = z 2 (8 + z 2 )K (z) + Az(A + z 2 )K x {z) . (11) 

We plot in Fig. 3 the mass spectrum for the same {tq,To} initial condi- 
tions as in the previous section. We have also shown the Drell-Yan spectrum 
for collision of two lead nuclei at yfs = 200A GeV, with structure functions 
obtained from the set I of Duke and Owens [|i~5]| . It is clear from this figure 
that our thermal QCD signal has little or no chance of shining through the 
hadron gas and Drell-Yan "backgrounds". Then, for a somewhat extreme 
viewpoint, we have also plotted the invariant mass spectrum for dileptons, 
taking initial temperature To as 500 MeV and the initial time tq = l/3To, 
appropriate for energies reached at RHIC for collisions involving two gold 



nuclei [IB]. Looking at Fig. 4, we immediately realize that the invariant 
mass window for seeing the emissions from the QGP are pushed to beyond 
2 GeV, once the full spectrum of dilepton yielding reactions is included. We 
note the discomforting feature that, even at such temperatures, the emissions 
from the QGP are neither much larger than the Drell-Yan contribution, nor 
are they significantly larger than the emissions from the hadronic matter. 
This demands a very careful analysis of the dilepton spectra before it can be 
used as a signature of QGP at RHIC energies. We do not plot the M T scaling 
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features of our solution with this high initial temperature: there, transverse 
expansion effects could certainly not be ignored. The transverse expansion 
influence on the net invariant mass spectrum also needs to be computed, as 
it will shorten the lifetime of the hadronic and mixed phases [jnj . 



3 Discussion 

We have addressed in this paper the issue of the proposed My scaling of 
the dilepton spectrum in high energy heavy ion collisions. We have found 
that the integrated lepton pair distributions would not keep any memory of 
the Mt scaling properties of the original plasma phase. This is because of 
the sheer size of the hadronic gas signal, and because the scaling hypothesis 
is violated strongly in the nonperturbative sector, owing to the structure of 
the timelike electromagnetic form factors. Of course, any calculation such as 
ours is an approximate one. However, we believe that the approximations 
we have made would have made the scaling easier to observe. 

Firstly, we have not attempted a quantitative calculation of transverse 
flow stemming from hydrodynamic effects. To reiterate, any flow features 
will spoil the scaling properties. If the quark-gluon signal is helped by dialing 
a high initial temperature in our calculations, this indicates (i) a somewhat 
longer lived QCD phase where flow will develop on its own and, (ii) the 
importance of the relaxation mechanisms. This brings us to our second point. 
In our transverse mass spectra we have neglected the Drell-Yan signal. It 
is easy to verify using simple expressions that the Drell-Yan signal violates 



13 



Mt scaling. This is easily understood, as one of the premises for scaling is 
thermal equilibrium. However, the simple Drell-Yan calculations might not 
be so relevant for the high energy collisions of heavy nuclei if the original 
parton distribution is modified during a nucleus-nucleus collisions. Such 
effects could surely modify the Drell-Yan estimates. Perhaps an adequate 
attempt to consider parton rescattering effects is provided by the Parton 
Cascade Model (PCM) 0. In this picture, the partons start in a completely 
nonequilibrium state and they are subsequently driven towards relaxation 
by continuous processes. In the PCM, the Mt scaling features have been 
analyzed The result there is that the scaling hypothesis is also strongly 
violated, even at the parton level. This state of affairs may owe to several 
facts. Among others, they are the lack of thermal and chemical equilibration, 
and the existence of additional scales in the problem, provided by infrared 
cutoff parameters and higher order QCD corrections. Also, processes that do 
not scale may also contribute Combining the above considerations with 
our own results, we are forced to the somewhat pessimistic conclusion that 
there is little hope that the observed lepton pair spectrum could show any 
evidence of "genuine" scaling. Conversely, any flatness in the transverse mass 
spectrum, plotted against transverse momentum, could not be interpreted 
unambiguously as a relic of a deconfined QCD phase. 
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Figure captions 

Fig. 1 The effective form-factor for dilepton yielding processes in hadronic 
matter when the rates obtained in Ref. [7] are written as Eq. (2.4). The 
dashed-dotted, dashed and solid curves represent the effective form- 
factor for T= 150, 200 and 250 MeV, respectively. The tight-dotted 
curve is the form- factor for pionic annihilation from Biagini et al. [11]. 
The loose-dotted curve represents the quark-ant iquark "form factor" . 

Fig. 2 Transverse-mass spectra for dileptons, for contributions from quark- 
annihilation (horizontal dashed line), from the full set of mesons (dashed- 
dotted curve), and from the pion gas approximation (solid curve). The 
temperature is 250 MeV and tq — 1 fm/c. The transverse mass is fixed 
at My = 2.6 GeV [8]. 

Fig. 3 Dilepton mass spectra. We show the contributions from quark- 
matter (dashed line), pion annihilation (dotted curve), Drell-Yan (dashed- 
dotted curve) processes and from all the hadronic processes considered 
in Ref. [7] (solid curve). The temperature and formation time are the 
same as in the previous figure. 

Fig. 4 Same caption as in Fig. 3, but for T = 500 MeV and r = 1/3T . 
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